function[sim] = muscle_vector_sim( musclevector_matrix )
% position: which normal vector
% i: which subject vector
% num: sim

[rows,cols] = size(musclevector_matrix);

for i = 1:cols
    [num,position] = max(musclevector_matrix(:,i));
    musclevector_matrix(position,:) = zeros;
    sim(i).musclevector_sim = num;
    sim(i).normal_vector_num = position;
    sim(i).subject_vector_num = i;
end
end

